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Abstract 

In 2005, Berenstein and Vazquez determined an open spin chain Hamiltonian de- 
scribing the one-loop anomalous dimensions of determinant-like operators correspond- 
ing to open strings attached to y = maximal giant gravitons. We construct the 
transfer matrix (generating functional of conserved quantities) containing this Hamil- 
tonian, thereby directly proving its integrability. We find the eigenvalues of this transfer 
matrix and the corresponding Bethe equations, which we compare with proposed all- 
loop Bethe equations. We note that the Bethe ansatz solution has a certain "gauge" 
freedom, and is not completely unique. 
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1 Introduction 



The discovery and exploitation of integrability in planar AdS / CFT has already led to many 
remarkable results [1], and may even ultimately lead to the solution of planar A/" = 4 su- 
persymmetric Yang- Mills (SYM), which is widely regarded as the "harmonic oscillator" of 
4-dimensional gauge theories. Much of the focus has been on the problem of computing the 
anomalous dimensions of single-trace operators of A/" = 4 SYM, which can be mapped to the 
problem of determining the eigenvalues of certain integrable closed spin-chain Hamiltonians, 
as observed in the seminal work of Minahan and Zarembo [2]. ^ However, progress has also 
been made on the problem of computing the anomalous dimensions of determinant-like op- 
erators, which can be mapped to the problem of determining the eigenvalues of certain open 
spin-chain Hamiltonians [6]. By the AdS/CFT correspondence, these two types of operators 
correspond to states of closed strings and open strings attached to D-branes, respectively. 

The simplest and most-studied open string/ chain example is the so-called Y — Q maximal 
giant graviton brane [7]. The one-loop SO{Q) scalar sector open spin chain Hamiltonian was 
found by Berenstein and Vazquez [8]. They also determined the (one-loop) boundary S- 
matrix, and showed that it satisfies the boundary Yang-Baxter equation (BYBE) [9, 10, 
11]. A breakdown of integrability in the SU{2) subsector at two loops was suspected in 
[12, 13]. However, Hofman and Maldaccna [7] subsequently showed that integrability is in 
fact preserved at two loops; and, based on SU {1\2) symmetry, they proposed an all-loop 
boundary S-matrix that satisfies the BYBE. (See also [14].) The scalar factor of the all- 
loop boundary S-matrix was proposed in [15]. Corresponding all- loop Bethe equations were 
proposed by Galleas in [16]. The bulk worldsheet Yangian symmetry discovered in [17], 
suitably generalized to the case of boundary scattering [18, 19], was used to construct all- 
loop bound-state boundary S-matrices [18, 20]. Finite-size corrections have been considered 
in [21, 22]. The classical integrability of the corresponding string sigma model with boundary 
has been investigated in [23, 24]. 

Despite the availability of all-loop results, there are still some important unresolved 
problems at one loop (where results are generally most solid) : the transfer matrix (generating 
functional of conserved quantities) containing the Berenstein- Vazquez Hamiltonian has not 

been constructed, and the corresponding eigenvalues and Bethe equations have not been 
determined. For the corresponding one-loop 5*0(6) scalar sector closed spin chain, such 
results were already obtained in the original Minahan-Zarembo work [2]. 

The purpose of our paper is to fill this gap: namely, to construct the one-loop open- 
chain transfer matrix, and to determine its eigenvalues. In this way, we directly prove the 
integrability of the one-loop Hamiltonian [8], and test the all-loop Bethe equations [16]. The 
simpler case of the SU{3) subsector was treated in [25, 26]. 

The outline of this paper is as follows. In Section 2, we recall the open spin chain 
Hamiltonian found in [8], and construct the corresponding transfer matrix. In Section 3 
we use the analytical Bethe ansatz to determine the eigenvalues of this transfer matrix and 
the corresponding Bethe equations. Surprisingly, the Bethe ansatz solution has a certain 

^We have in mind here "long" operators. For operators of finite length, there are finite-size corrections 
[3, 4, 5]. 
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"gauge" freedom, and is not completely unique. One of the sets of one-loop Bethe equations 
that we find is consistent with the all-loop equations. We conclude in Section 4 with a brief 
discussion of our results. 



2 Construction of the transfer matrix 

In [8] Berenstein and Vazquez identified an open 5*0(6) spin chain that describes the one- loop 
anomalous dimensions of determinant-like operators corresponding to open strings attached 
to y = maximal giant gravitons. The space of states is ^ 

1 L L+1 

®---C^ . (1) 

That is, the vector spaces of the "bulk" sites (labeled 1, . . . , L) all have dimension 6, while 
the vector spaces of the left and right "boundary" sites (labeled and L -|- 1, respectively) 
have dimension 5. The Hamiltonian is given in Eq. (2.15), which we rewrite as 

L-l 

H = Qlh,,,Ql + Ql) + hi,i^^ + Ql+M,L+iQl+i + {1- Ql+i) , (2) 

1=1 

where is the bulk two-site Hamiltonian 

hi,i+i = -^JCi^i+i + - Vi,i+i . (3) 

We have relabeled the sites to run from to L -|- 1 (instead of 1 to L) ; we have relabeled Z 
and Z hj Y and Y, respectively; and for simplicity, we have set the coupling constant A = 1. 
We note that Q'^ is the projector 

Q*|0) = O, Q*|(^) = |(^)for(^^0. (4) 

In the standard basis \a) — Ca with a — 1, . . . , 6 (elementary 6-dimensional vectors with 
components [ca]- — Sa,i), the matrices V and /C are given by 

6 6 
7^ = Cab <8) eta , /C = Cafe (8) Cab , (5) 

a, 6=1 a,b=l 

where Cab is the standard elementary 6x6 matrix whose (a, b) matrix element is 1, and all 
others are zero; i.e., [eaf,] . . = SaiSbj. 



It is convenient to change to a new basis. 


|2) = 


-\W)- 




- ie2) , 




|3) = |z) = -i=(e3 + ze4). 


|4) = 










|5) = \Y)^^{e, + ie,), 


|6) = 






■ iee) . 


(6) 



^Following [7], we define the origin of tlie spin chain at site instead of site 1. 
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Let U be the unitary operator which implements this change of basis, 

\a) = U\a) . 

Its matrix elements in the original basis Uba — {b\U\a) — {b\a) are given by 



(7) 



U = 



V2 



/I 1 

i —i 



V 



1 1 

i —i 



1 1 
i —i / 



(8) 



where matrix elements that are zero are left empty. 

In this new basis, V does not change, but JC does change: 

p = {U^ ®U^)V{U®U)=V, 

iC = {U^ ®U^)}C{U^U). (9) 
Moreover, in this new basis, the projectors and are given by the diagonal matrices 

= diag(l,l,l,l,0,l) = l-|F)(F|, 

= diag(l,l,l,l,l,0) = l-|F)(F|. (10) 
We therefore arrive at the following explicit matrix representation of the Hamiltonian 

L-l 



(11) 



+ — Vi,i+i 



1=1 

where the bulk two-site Hamiltonian is given by 

hi,i+i 

and the boundary terms are given by 

Hbt = Qoho^iQo + {I - Qo) > 

= Qj^^lhL,L+lQl+i + {I - Ql+i) ■ 



The boundary terms have the property 

^bt — Qo ^bt Qo ' — Q\+1 Q\+1 ■ 



(12) 

(13) 
(14) 

(15) 



We drop the null rows and columns in these matrices that are due to the projectors. 
Hence, i?^ and should be understood as 30 x 30 matrices acting on x and x C^, 
respectively, as indicated in (1). 
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We now address the problem of constructing the corresponding transfer matrix that 
contains this Hamiltonian as well as the higher conserved charges. According to Sklyanin 
[10], in order to construct an open-chain transfer matrix, we need an i?- matrix that gives 
the bulk two-site Hamiltonian; and also left and right i^-matrices, which give the left and 
right boundary terms, respectively. ^ 

2.1 R(u) matrix 

We recall that the Yang-Baxter equation is given by 

i?12(Ml - U2) RlsiUl) ^23(^2) = ^23(^2) ^13(^1) i?12(Ml - U2) ■ (16) 

An S'0(6)-invariant solution which acts on is given (in the new basis) by [27, 2] 

R(u) = \u{2u + 2 - n)/ - (2m + 2 - n)V + 2u)c\ , (17) 

n — 2 I J 

with n — 6. This i?-matrix indeed produces the bulk two-site Hamiltonian (12), since 

hi,i+i^Vi,i+iR[i^M + lhi+i: (18) 
where the prime denotes differentiation with respect to the spectral parameter u. 

2.2 K~{u) matrix 

The right matrix should give the right boundary term (14) in the Hamiltonian. The 
K-matrix must therefore be operator-valued, rather than c-number valued. Indeed, the 
K-matrix K~{u) (which acts on C^) must satisfy the right "operator" BYBE 

^12(^1 - U2) K{^{ui) ^12(^1 + U2) ^^23(^2) 

= K~^{u2) Ri2{ui + U2) K~^{ui) Ri2{ui - U2) , (19) 

where the i?- matrix is given by (17). 

We claim that the needed X-matrix is given by 

K{2{u) = QlRi2{u)Ki{u)Ri2{u)Ql , (20) 

where K{u) is the 6x6 diagonal X-matrix [28, 29] 

K{u) = diag(l — u,l — u,l — u,l — u,u — l,u + 1) , (21) 



^In order to avoid confusion, it may be worth noting that these R and K matrices, even though they 
satisfy bulk and boundary Yang-Baxter equations, have no direct relation to the bulk and boundary S- 
matrices discussed in the Introduction. This fact is evident in the more familiar case of the ferromagnetic 
spin-1/2 XXX Heisenberg chain: while the bulk S-matrix is a U{1) phase, the i?-matrix is an S'f/(2)-invariant 
4x4 matrix. 
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which satisfies the standard (not operator) BYBE 

RuiUi - U2) Ki{Ui) Ri2{Ui + U2) k2{u2) 

= k2{u2) Rl2{Ul + U2) ki{ui) Ri2{Ui - U2) . (22) 

We note that operator /T-matrices of the general type (20) were introduced by Frahm and 
Slavnov [30], who called them "projected" /^-matrices. The fact that i^-matrices of this 
type are needed to construct the SU{'i) subscctor of the F = spin chain was noted in [26]. 

Indeed, we have verified that (20) does satisfy the operator BYBE (19), and produces 
the right boundary term (14), 

H^t^\K-\Q)-r2I. (23) 

We also note that 

X-(0) = 7. (24) 



2.3 K+{u) matrix 

The left fC-matrix should give the left boundary term (13) in the Hamiltonian. For the 
SUi^i) case [25] and in fact for the general SU{N) case, we found in [26] that the needed 
left X-matrix can be obtained from the right X-matrix using Sklyanin's [10] "less obvious" 
isomorphism ^ 

K+{u) = kiu)tT2Vi2Ri2{-2u-7j)K^^{u), (25) 
where rj appears in the crossing-unitarity relation 

Ri2iuY'Ri2{-u-vf' ^I, (26) 

and k{u) is an arbitrary scalar factor. 

Fortunately, the same trick works also for the 5*0(6) case. Indeed, the 5*0(6) i?-matrix 
(17) satisfies the crossing-unitarity property (26) with r] = —4. We have verified that the 
matrix K'^{u) given by the isomorphism (25), which acts on (7^ (g) (7^, satisfies the left 
operator BYBE 

Rui-ui + U2) K+iuiY' Rui-ui -U2-v) K+{u2y' 

= KUU2Y' Ri2{-ui -U2-v) KUulf' Rui-u, + U2) . (27) 
And, most importantly, this X-matrix produces the left boundary term (13), 

I k i^+o(0) «i(0)P„i + VaiR'AO)) + tra K+ '(0)1 + ll , (28) 
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''While in [25] the left if- matrix is defined to act on VqVi, here we define the left if- matrix to act on ViVq; 
i.e., the two if- matrices are related simply by permutation of the vector spaces Vo and Vi. The isomorphism 
(25) evidently gives directly the latter form. 
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where the trace is over a 6-dimensional auxihary space, which is discussed further below. 
We have fixed the scalar factor in (25) to be k{u) = [u{u — 2)(2u — 1)]^^ in order to cancel 
corresponding terms that appear in the numerator. Finally, we note that 

tr„X+(0) = -3/. (29) 



2.4 Transfer matrix 

Having gathered all the necessary ingredients, we are now ready to assemble them to form 
the transfer matrix. We introduce a 6-dimensional auxiliary, denoted by a, and consider 
operators on the enlarged vector space (cf. (1)) 

a 1 L L+1 

®---C^ . (30) 

We define the monodromy matrices T and T by 

Tai...L{u) = Ral{u) ■ ■ ■ RaL{u) , fai...L{u) = RaL{u)---Ral{u), (31) 

where R{u) is given by (17). The transfer matrix is given by [10] 

t{u) = tr„ K^^iu) T,,...l{u) K-^^,{u) t^.-.M , (32) 

where K^{u) and K^{u) are given by (20) and (25), respectively. 

By construction, the transfer matrix has the fundamental commutativity property 

[tiu),t{v)]^0, (33) 

and contains the Hamiltonian (11), 



(34) 



The relations (33) and (34) directly imply the integrability of the Hamiltonian. Higher 
conserved quantities can be obtained from higher derivatives of the transfer matrix at u = 0. 

We note that the transfer matrix t{u) is crossing invariant up to a scalar factor, 

u 

- = ^^(«) ■ (35) 
Equivalently, defining a rescaled transfer matrix t{u) by 

i{u) = ut{u), (36) 
we see that this rescaled transfer matrix is exactly crossing invariant, 

t{2 -u)^ t{u) . (37) 
The Hamiltonian is evidently related to the second derivative of i{u) at w = 0. 
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3 Analytical Bet he ansatz 



The commutativity property (33) implies that it is possible to find eigenstates |A) of the 
transfer matrix t{u) (32) which are independent of u, 



t{u) |A) = A{u) |A) . 



(38) 



We turn now to the problem of determining the eigenvalues A{u). We proceed by the 
analytical Bethe ansatz approach, along the fines in [2, 25]. We choose as a reference state 



|0) = 



/1\ 








/1\ 








/1\ 







(39) 



which is a ground state of the Hamiltonian (11) with eigenvalue 0. We denote the corre- 
sponding eigenvalue of the transfer matrix by Ao{u), 



tiu) |0) = Ao(«) |0) . 



(40) 



On the basis of results for L = 0, 1,2, we obtain the following conjecture for the vacuum 
eigenvalue 

1 



a(K)(K-2)2^(H-l)'^ + 6(K)(H-l)'^K'^ + 4c(ii)(K-2)^^K^^ , (41) 



\2L 



2L„,2L 



Ao{u) 
where 

a{u) = {u-2)\u-lf{2u-3)\ 

b{u) = {u - 2){u - l)^u%2u - If , 

c{u) = 4(ii - 2) V - 1) V , 

d{u) = 16(2m-3)(2m-1). 

A general eigenvalue should be a "dressed" vacuum eigenvalue, ^ 



\2L„.2L 



(42) 



ii-l)^V^ 



Qi(^-i) 

Qi{u-I) 



+ (m-2)^V^ ci(«) 



Qi{u - l)Q2{u - l)Q-i{u - 1) 



+ C2[U 



Qi{u-^)Q2{u-2)Qs(u-2) 
'Qi(u-1)Q2{u-1)Qs{u-1) 



Q2{u)Q3{u-2) Q2iu-2)Qs{u 

-r C3[U)-——. TvFTl TT "T ''4l,"J 



Q2{U-1)Q3{U-1) 



Q2{u-l)Qs{u-l) 



(43) 



^For the corresponding closed-chain result, see (4.28) in [2] and references therein. 
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where 

ma 

Qa{u) = Yl{u-iUa,k){u + iUa,k) , = 1,2,3, (44) 

k=l 

which have the property Qa{—u) — Qaiu)- The functions ci{u), . . . , C4{u) must satisfy 

ci{u) + C2{u) + C3{u) + C4{u) = Ac{u) , (45) 

but are otherwise still to be determined. 

The eigenvalues of the Hamiltonian (11) now follow from (34) and (43), 

1 1 7 1 

which is the same expression that one obtains for the corresponding closed chain [2]. 
The crossing symmetry (35) implies a corresponding property of the eigenvalues 

7/ 

A(2 -u)^ 7— > (47) 
which in turn imphes the constraints 

ci(2 -u) = :r—C2{u) , C3(2 -u)^ 7r^C4{u) . (48) 

Zl Hi Zl Hi 

The Bethe equations for the zeros Ui^k of the function Qi{u) (44) follow from the fact 
that A(k) in (43) is analytic aX, u — iui^k + which implies 

wi,fc + |y^_ a(mi,fe + i) (5i(wi,fc + 1) Q2(?wi,fc - |) Qsl^Mi^fc - |) 



iui,k -\j ci(mi,fc + \) Qi{iui^k - 1) Q2{iui,k + |) Q^{iui,k + \) 

We should obtain the same Bcthe equations by considering instead the poles aX,u — mi,A; + 1, 
which implies the constraint 

a{ul C2{u^_ (50) 

Similarly, by considering the poles at m = iu2,k + 1, we obtain the Bethe equations 

^ _ Ci(m2,fc + 1) Q2iiU2,k + 1) Qi{iu2.k - I) . s 

C4{iU2,k + 1) Q2{iU2,k - 1) Ql{iU2,k + \) ' 



and the constraint 

Cl{u) ^ C3{u) 
C^iu) C2{u) 



(52) 
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Finally, by considering the poles at u — iu^^k + 1, we obtain the Bethe equations 



i) 



C3iiu3,k + 1) Qsiius^k - 1) Qi{ius,k + o) 



(53) 



and the same constraint (52). It can be shown that the requirement that the eigenvalues 



(43) be analytic at u 



1 3 
2' 2 



(where d{u) has zeros) does not lead to further constraints. ^ 



The constraints (45), (48), (50), (52) do not uniquely determine the functions Ci{u), . . . , C4{u). 
Moreover, the problem of determining the transfer-matrix eigenvalues A[u) in the T — Q form 
(43) does not have a unique solution. Indeed, consider the following ansatz for these functions 



ci{u) 

C3{u) 
C4{u) 



m"(m-1)^-"(m-2)^(2m-3)\ 



u 



3+m/ 



{u - ly-^'iu - 2y+"{2u - l){2u - 3) , 



{u - If-'^iu - 2)^+"'{2u - l){2u - 3) , 



(54) 



where n and m are integers. All the constraints (45), (48), (50), (52) are then satisfied for 
the following four sets of (n, m) values 



case I : 


n 


= 3, 


m 


= 0, 


case Ila : 


n 


= 5, 


m 


= 0, 


case lib : 


n 


= 5, 


m 


= 2, 


case III : 


n 


= 7, 


m 


= 2, 



(55) 



which we have designated clS CclSGS 1, Ila, lib, III, respectively. 

The Bethe equations (49), (51), (53) can now be rewritten in the more familiar form 

mi m2 

Ue2(Ml,fc - Mlj) e2(Ml,fe + Mij) JJe_i(Mi,fe - Msj) e_i(Mi,fe + M2j) 



ei(Mi,fc) 



2L+n-l 



X e_i(Mi,fc - M3j) e_i(Mi,fc + M3j) , A; = 1 , . . . , mi , 



(56) 



m2 



e2{u2,kT H e2(M2,fc - U2,j) e2{u2,k + U2,j) 



X JJe_i(ii2,ifc-iiij)e_i(ii2,ifc + iiij), A; = l,...,m2, (57) 



^It may be possible to derive further constraints from the requirement that the eigenvalues of the transfer 
matrix obtained by fusion in the auxiliary space also be analytic. However, such calculations would be 
difficult, and we shall not pursue them here. 
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1 



e2(M3,fc) 



n—m—3 



X 



Yl e-i{u3,k - uij) e_i(u3,jk + uij) , /c = 1 , . . . , ma , 



(58) 



where 



en{u) 



u + f 



(59) 



For cases I and III, the Bethe equations are symmetric under the interchange of Bethe roots 
of types 2 and 3 (i.e., U2,i -H- ms,^), as in the closed-chain case [2]. The Bethe equations for 
cases Ila and lib transform into each other under the interchange of Bethe roots of types 2 
and 3. 

At least for small values of L, the Bethe ansatz solutions corresponding to each of the 
cases (55) are complete; i.e., they reproduce all the transfer-matrix eigenvalues, of which 
there are 6'^5^ according to (1). Indeed, following the approach described in Appendix B of 
[25], we have verified the completeness for L = and L = 1, for which the total number of 
states is 25 and 150, respectively. The results for case I are summarized in Table 1 and Table 
2, respectively. For cases II and III, the Bethe roots of type 1 (i.e., {ui^k}) describing a given 
eigenvalue are the same as for case I, but the Bethe roots of types 2 and 3 are different.^ This 
does not lead to differences in the energy (46) or higher conserved quantities, which depend 
only on the Bethe roots of type 1. Hence, the various cases can be regarded as equivalent. 
Gauge-like transformations relating these cases are discussed in Section A. 

For case I, the above one-loop Bethe equations coincide with those that we previously 
conjectured, and obtained from the all-loop Bethe equations [16] by performing the weak- 
coupling limit and then reducing to the S0{6) sector, which are given in Eqs. (3.1) and 
(4.20) in [25], respectively. For cases II and III, the Bethe equations are evidently slightly 
different: in (56) the term ei{ui^k) has a different power; and (57) and/or (58) contain an 
additional factor el- 

4 Discussion 

We have shown that the Berenstein- Vazquez Hamiltonian (11) is contained in the commut- 
ing transfer matrix (32), which directly implies the integrability of this Hamiltonian. One 
key ingredient is the left matrix (20), which is of the projected type [30]. Another key 
ingredient is the right X-matrix (25), obtained through a seldom-used isomorphism noted in 
[10]. We have also found expressions (43), (44), (54), (55) for the eigenvalues of this transfer 
matrix in terms of roots of the Bethe equations (56)- (58). We have therefore completed the 
generalization of important closed-chain results of Minahan and Zarembo [2] to the Y = 
open chain. 

^Generally, there are fewer Bethe roots of types 2 and 3 for cases II and III in comparison with case I. 
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We have seen that the Bethe equations based on the vacuum eigenvalue (41) are not 
uniquely fixed - a certain "gauge" freedom exists. A similar phenomenon may occur for other 
intcgrable spin chains with symmetry group of rank greater than one. (This phenomenon is 
distinct from the well-known duality transformations of supersymmetric spin chains, see e.g. 
[31].) For case I in (55), the one- loop Bethe equations (56)- (58) agree with those obtained in 
[25] from the all-loop Bethe equations [16] by performing the weak-coupling limit and then 
reducing to the SO (6) sector. Hence, our results provide support for those all-loop Bethe 
equations. 

The analytical Bethe ansatz approach that we followed here is heuristic, and requires 
making some assumptions, in particular (54). It would be useful to carry out a more rigorous 
analysis based on nested algebraic Bethe ansatz. For the 0{2N) closed spin chain, a nested 
algebraic Bethe ansatz was developed in [32] . One would need to generahze that approach to 
the open chain with projected X-matrices considered here. Besides the difficulty of working 
out the necessary commutation relations, one can foresee the following further difficulty: after 
the first level of nesting, the "bulk" quantum space will have 4 dimensions (two less than the 
original 6 dimensions), while the "boundary" quantum spaces will still have 5 dimensions. 
That is, the bulk and boundary quantum spaces will not have the same dimension. The 
reduced transfer matrix will therefore not be given by the tensor product of two usual 
transfer matrices, as happens at the final stage in the GL{N) case [26]. 

We have restricted our attention here to the Y = Q maximal giant graviton brane, which 
is the simplest known integrable open string/chain example. It may be interesting to consider 
the corresponding problem in more complicated examples, such as the Z = maximal giant 
graviton brane [7]. 
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A "Gauge" transformations among the cases 

Let us denote by cJ(m), cf^{u), etc. the functions (54) for cases 1, 11a, etc. in (55), respectively. 
It is easy to see that the functions for case Ila are related to the corresponding ones for case 
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I as follows 



u 



u-l 



u-2 
u-l 

u-2 
u-l 



u 



c\{u 



dyiu 



c\{u 



c\{u 



u-1^ 

Similarly, the functions for case lib are related to the corresponding ones for case I as follows 

2 



Jib/ 



U 



u-l 



u-2 
u-l 



c\{u 



C2[U 



u 



« \ 1/ 



u-2 
u-l 



c\(u 



Finally, the functions for case III are related to the corresponding ones for case I as follows 



u 



c'^\u) 



u-l 



u-2 



c\{u 



c^{u 



Co m 



u-l 

u''{u-2 f I 
{u-l) 



(60) 



(61) 



(62) 



Let us now consider the expression (43) for A(-u), and identify the functions q(-u) there 
with <^^{u). Our observation is that, by making in that expression the simple transformation 



Q^{u) u^Qz{u) , 



(63) 



we obtain, in view of (60), the corresponding expression for A(m) in terms of cf^{u). Similarly, 
by making instead the transformation 



Q2{u) u^Q2{u) , 



(64) 
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we obtain, in view of (61), the expression for A(m) in terms of cY^{u). Finally, in view of 
(62), by making both transformations (63) and (64), we obtain the expression for A{u) in 
terms of c™(m). 

Similarly, let us consider the Bethe equations in their original form (49), (51), (53), and 
identify the functions Ci{u) there with cl{u). By making the transformations (63) and/or 
(64), we obtain the Bethe equations for the other cases. 

In short, (63) and (64) are sorts of (discrete) "gauge" transformations that relate the 
four cases (55) of Bethe ansatz solutions. In our numerical investigations, we have observed 
that for case I there arc generally more zero Bethe roots of types 2 and 3 in comparison 
with cases II and III, which is consistent with the presence of additional factors of in the 
corresponding Q functions. 

Related transformations were briefly discussed in [26]. There, the transformations in- 
volve also Qi{u), and therefore, relate Bethe ansatz solutions based on reference states with 
different energies. 



deg 








9 
6 
4 
4 
2 


1/2 
VS/6 
V3/2 
±i/2 














Table 1: Degeneracy and Bethe roots for case I with L — 0. 
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deg 




{'^2,k} 


{U3,k} 


16 


— 


— 


— 


16 







— 


16 


\/3/2 





— 


9 


1/2 


— 


— 


9 


1/2 








9 


(v^+l)/2 








9 


(a/2-1)/2 








6 


±i/2 





- 


10 


±i/2 








14 


0.230955, 0.668326 








14 


0.716015 ±0.512521i 








14 


VS/6,VS/2 


v^/3 





4 


0.41551L 1.15211 


1 


1 


1 


1/2, ±V2 


0,i 


0,i 


1 


1/2, 0.479716 ± 0.971633i 


0, 1.38848i 


0, 1.38848i 


1 


1/2, 0.208963, 1.02227 


0, 0.767271 


0, 0.767271 


1 


1/2, 0.414496 ±0.502211i 


0, 0.812907i 


0, 0.812907i 



Table 2: Degeneracy and Bethe roots for case I with L — 1. 
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